Effective numerical methods for simulating diffusion on a spherical surface in three dimensions

In order to construct an algorithm for homogeneous diffusive motion that lives on a sphere, we consider the equivalent process of randomly rotating spin vector constant length. By introducing appropriate sets random variables based cross products, families methods with increasing efficacy exactly preserve modulus every realisation. This is done by exponentiating antisymmetric matrix whose entries are these Gaussian in simplest case.